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Abstract 

In this paper I explore whether a vector field can be the origin of the present stage of cosmic 
acceleration. In order to avoid violations of isotropy, the vector has be part of a "cosmic triad" , 
that is, a set of three identical vectors pointing in mutually orthogonal spatial directions. A triad 
is indeed able to drive a stage of late accelerated expansion in the universe, and there exist tracking 
attractors that render cosmic evolution insensitive to initial conditions. However, as in most other 
models, the onset of cosmic acceleration is determined by a parameter that has to be tuned to 
reproduce current observations. The triad equation of state can be sufficiently close to minus one 
today, and for tachyonic models it might be even less than that. I briefly analyze linear cosmological 
perturbation theory in the presence of a triad. It turns out that the existence of non-vanishing 
spatial vectors invalidates the decomposition theorem, i.e. scalar, vector and tensor perturbations 
do not decouple from each other. In a simplified case it is possible to analytically study the 
stability of the triad along the different cosmological attractors. The triad is classically stable 
during inflation, radiation and matter domination, but it is unstable during (late-time) cosmic 
acceleration. I argue that this instability is not likely to have a significant impact at present. 
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I. INTRODUCTION 

A combination of different cosmic probes that primarily involves supernova data 
suggests that the universe is presently undergoing a stage of accelerated expansion. Little 
is known about the origin of this stage of cosmic acceleration Q]. It might be related to a 
breakdown of general relativity on large scales , or it can be the effect of "dark energy" 
(|, a negative pressure component that causes the universe expansion to accelerate. The 
simplest possibility is that dark energy merely is a (tiny) cosmological constant. If dark 
energy is dynamical, it is mostly assumed to be a scalar field, quintessence Q. 

On purely phenomenological grounds, within the context of general relativity, dark energy 
might be characterized by its equation of state wng, its speed of sound c 2 s and its anisotropic 
stresses Q]. Conventional quintessence models [J have an equation of state — 1 < wde, a 
speed of sound c 2 = 1 and no anisotropic stresses, whereas in k-essence models P, [kJ the 
speed of sound is in general different from one. The range of possible phenomenologica 
properties of dark energy is not exhausted by the former models. Phantom dark energy 111 | 
has an equation of state Wde < — 1, and thus violates the dominant energy condition [12] . 
The origin of such an equation of state is the "wrong" sign kinetic term of the phantom scalar 
field. Because of that, phantom dark energy is quantum- mechanically unstable upon decay 
of the vacuum into (positive energy) gravitons and (negative energy) phantom particles 
12, 13. Hence, if future preciser observations confirm the current trend and favor a dark 

fin 

energy equation of state smaller than minus one |l|, Il4j . alternative (viable) models are 
needed to account for such a value 1 jl(| . 

In this paper I explore whether dark energy can be vector-like. Vector-like dark energy 
turns out to display a series of properties that make it particularly interesting phenomeno- 
logically. On one side, it can also violate the dominant energy condition, w^e < — 1> while 
possessing a conventional kinetic term. On the other side, it has non-anisotropic stress per- 
turbations and it leads to violations of the decomposition theorem Q], i.e. the decoupling 
of scalar, vector and tensor cosmological perturbations. As we shall see, other interesting 
features include the existence of tracking attractors, that render the evolution of a vector 
field in an expanding universe rather insensitive to initial conditions. In spite of this attrac- 



1 The vacuum-driven metamorphosis model of Parker & Raval (lfij seems to be experimentally ruled out at 
the 99.5% confidence level [J. 



2 



tive feature, the onset of cosmic acceleration is determined by a parameter in the Lagrangian 
that has to be properly adjusted, as in most other models. In that respect, these forms of 
vector dark energy are similar to quintessence trackers 1(3, Q] . 

Non-gravitational interactions are known to be mediated by vector fields. In addition, 
from a four-dimensional point of view, certain components of higher-dimensional metrics 
behave like vectors [3]. It is hence natural to study the evolution of vector fields in a 
cosmological setting. However, the existence of a spatially non-vanishing vector breaks the 
isotropy of a Friedmann-Robertson- Walker (FRW) universe. From the point of view of 
gravity, such a breaking manifests itself in anisotropic stresses caused by the vector. If dark 
energy is such a vector, as long as it remains subdominant, this violation is likely to be 
observationally irrelevant j^j]. Once dark energy comes to dominate though, one would 
expect an anisotropic expansion of the universe, in conflict with the significant isotropy of 
the cosmic microwave background (CMB) 2l|. Because of that, in this paper I consider a 
"cosmic triad", i.e. a set of three equal length vectors that point in three mutually orthog- 
onal spatial directions. Remarkably, the existence of a triad turns to be compatible with 
spatial isotropy, at least from the point of view of gravity. While the triad guarantees the 
isotropy of the background, it does not automatically imply the isotropy of its perturbations. 
Eventually, it might be even necessary to introduce fields that explicitly violate rotational 
symmetry as there appear to be hints of (statistical) anisotropy in the CMB fluctuations 
[22| . Along these lines, I speculate below that a triad could provide a link between cosmic 



acceleration and some of the anomalies observed in the CMB radiation 
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23]. 



Mainly because they single out spatial directions, vector fields have received compara- 
tively little attention in cosmology. Ford has proposed an inflationary model where a vector 
is responsible for a stage of inflation 24]. Our treatment here is to some extent similar 
to his proposal. Jacobson and Mattingly have studied the dynamics of a vector with of 
Sxed length, with the specffic purpose of studying violations of Lorent. invariance Q. A 
vector-like form of quintessence has been also considered by Kiselev [26], though his vectors 
significantly differ from ours. Zimdahl et al. have suggested that a (timelike) vector force 
could be responsible for the present acceleration of the universe J27 1. Also, it has been noted 
that the addition of higher-order powers of the Maxwell field-strengt 
of an electromagnetic field might cause the universe to accelerate 



1 to the Lagrangian 
The literature on 

magnetic lielxls in the early universe is more extensive, see [29] and references therein. 
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II. VECTOR DARK ENERGY 



Consider a set of three self-interacting vector fields A" Strictly speaking, this is really a 
set of three one-forms, but I shall call them vectors. Latin indices label the different fields 
(a, b, . . . = 1 . . . 3) and greek indices their different spacetime components (/i, v, . . . = ... 3). 
As we shall see below, this number of vector fields is dictated by the number of spatial 
dimensions and the requirement of isotropy. We would like to study the dynamics of such a 
"cosmic triad" in the presence of gravity. Consider hence those vectors minimally coupled 
to general relativity, 

3 



S = d 4 Xy/^g 



R 
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where F£ v = d^A a v - d u A^ and A a2 = g^A^A". The action (0) thus contains three identical 
copies of the Lagrangian of a single vector field. The term V(A a2 ) is a self-interaction 
that explicitly violates gauge invariance. For completeness, in the Appendix I show how a 
triad could naturally appear from a gauge theory with a single SU(2) gauge group. In the 
following, I use the Einstein summation convention throughout, where a sum is implied only 
over indices in opposite positions. The indices that label the different vectors are raised and 
lowered with the flat "metric" 5 a b- 

The kinetic term F 2 in the action (0) is not unique in the following sense. Up to bound- 
ary terms, V ' V A^J V A^ is the only additional diffeomorphism invariant quadratic term that 



contains two derivatives of the vector A^ 
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301 ] . Because the dynamics of vectors known to 
occur in nature are described by a Maxwell term, I consider a F 2 term only. Additional cou- 
plings of the vector are constrained by tests of gravity j^] and limits on possible violations 
of Lorentz symmetry 2 [3] . To avoid such violations, I assume that the matter Lagrangian 
C m only depends on the metric S „„ and on the remaining matter fle.ds but not on the 
cosmic triad A a . In that respect, the triad is analogous to conventional quintessence |34J. 
Varying the action (0) with respect to the metric one obtains the Einstein equations 



2 Note that in a FRW universe Lorentz symmetry is (spontaneously) broken anyway, in the sense that 
there are non-vanishing vector fields, like the gradient of the Ricci scalar or the CMB temperature, that 
define a preferred direction. Such a breaking could be detected by non-gravitational experiments if the 
non- vanishing vector directly couples to matter See also [ij] for a clear discussion of the relation 

between coordinate invariance, Lorentz invariance and isotropy. 
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G/J.V — SttGT^, where the energy momentum tensor of the triad is given by 



(2) 



This energy momentum tensor is the sum of the three different energy momentum tensors 
of the decoupled vectors, ^'T^ = ^ a ^T^, neither of which is of perfect-fluid form. By 
varying the action with respect to the vectors A", one obtains their equations of motion, 

d,(V^F a n = 2y/=g-^A w . (3) 

The four-divergence of the last equation yields a constrain equation for the vector. As a 
consequence, each vector has three dynamical degrees of freedom, as it should. 

We shall study the dynamics of these vectors in a flat, homogeneous and isotropic FRW 
universe with metric 

d s 2 = -dt 2 + a 2 (t)dx 2 . (4) 

An ansatz for the vectors that turns to be compatible with the symmetries of this metric 
(homogeneity and isotropy) is 

A\ = 5\A{t)-a. (5) 

Hence the three vectors point in three mutually orthogonal spatial directions and they share 
the same time-dependent length, A a2 = A a ^A a ^ = A 2 (t). Substituting the ansatz and 
the metric (jlj) into the vector equations of motion Q I find 

A + 3 HA + ^ + ^A + d ^ , (6) 

where a dot means a derivative with respect to cosmic time t. Note that the 0-component 
of equation (jH} forces Aq to vanish, as in the ansatz (JHJ)- Substituting the metric (jlj) into 
the Einstein equations one obtains 

H 2 = ^P, (7a) 

^ = -^(p + 3 P ), (7b) 
a S 

where H = a/ a is the Hubble "constant" . The energy density of the universe is p = — T° and 
its pressure is defined by p-8j = T l j, where i and j run over the spatial spacetime components. 
Note that the energy momentum tensor has to be compatible with the symmetries of the 
metric. For the FRW metric (J3J), G°i = 0, so that T°j should also vanish. It can be easily 



verified that this is indeed the case for the ansatz The requirement of isotropy is 

non-trivial for a single vector, since its energy momentum tensor is 

(a) r, = Q(i + HA) 2 - Vj ^ + (^^A 2 - (A + HA) 2 ^ ^ • 6 a r (8) 

Although this energy momentum tensor is diagonal, its value along the j = a direction is 
different from the one along the directions perpendicular to it. Nevertheless, the total energy 
momentum tensor ^T^ u = ^ a has isotropic stresses, and the corresponding energy 

density pa and (isotropic) pressure pa are given by 

PA = \(A + HA) 2 + W(A 2 ), (9a) 

PA = \{A + HA) 2 - 3V(A 2 ) + 2^A 2 . (9b) 

Note that the equation of motion (JUJ) can be also derived from the condition 
Pa + 3H{ Pa +Pa) =0. 

To conclude this section let me point out a remarkable property of the cosmic triad. 
Namely, its equation of state w a = Pa/ Pa can become less than —1 (with a positive energy 
density) if dV/dA 2 is negative 35j] . Because a mass term for a vector A^ has the form 
V{A 2 ) = '^j-A^A^ 1 , I call such vectors "tachyonic". Tachyons (particles of negative squared 
mass) are usually associated with instabilities. In many cases, an instability merely signals 
the tendency of the system to evolve. In cosmology, those instabilities are not particularly 
terrible. In fact, a universe in stable equilibrium would be pretty lame, as it would not 
even expand. In the absence of gravitational instability structures would not form, and 



without the (effective) tachyonic mass of a scalar, it wou 



invariant spectrum of perturbations during inflation 
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oe quite difficult to seed a scale 
Scalar tachyons 3 indeed have 



been widely considered in the literature. Other forms of instability are more worrisome, like 
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3- 



By 



the quantum-mechanical instability of the vacuum in the presence of a phantom 
simple analogy, any form of instability associated with a tachyonic vector is not expected 
to be of this second kind, as the vectors have a conventional kinetic term. This question is 
not only of academic interest, since analyses of observational data (marginally) favor a dark 
energy equation of state w^e < — 1 [llQ 

I shall not deal here with the quantum mechanics of tachyonic particles, which even 
for scalars is not free of problems. Nevertheless, I also want to present some arguments 



By a "scalar tachyon" I mean a scalar field ip with a convex potential, d 2 V/dip 2 < 0. 
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that suggest that a tachyonic vector might be similar to a phantom field 38|. One of the 
arguments goes back to the Stueckelberg theory of massive vectors j^], 4o|. Consider for 
simplicity a vector field in Minkowski space, 

L = - K ^-A»A». (10) 

The field is massive for k — 1 and tachyonic for k = —1. Upon the substitution 
— > + m^dftS, the Lagrangian (fTUj) reads 

L = -F^F - ^-A^ - KmA^S - fasWS. (11) 



Note that (fTTj) contains an additional scalar, the Stueckelberg field S. For a massive vector 
(k = 1) S has a conventional kinetic term, but for a tachyonic vector (k = —1), its kinetic 
term has the "wrong" sign, like the one of a phantom |ll|. However, the additional field S 
turns to be constrained in the quantum theory. Even though it describes a massive vector, 
the Lagrangian (JTTJ) has a gauge symmetry, A^ — > A^ + <9 M A, S — > S 1 — mA, where A is a 
scalar function. Upon quantization, this gauge freedom is fixed by imposing the Stueckel berg 
subsidiary condition d^A^ — mS = 0, which relates S and the divergence of the vector [40]. 
Thus, strictly speaking, the field S is not simply a phantom scalar in the conventional 
sense. On the other hand, there are other properties that suggest phantom-like behavior of 
a tachyonic vector, like the opposite sign of the propagator in the high-momentum limit, or 
the opposite sign in front of the squared longitudinal momentum in the Hamiltonian |4l| . 

In the time-dependent situation we are dealing with, where the triad vectors have a 
non-vanishing expectation value, the issue is slightly more complicated. Consider quantum 
fluctuations 5A^ around one of the triad vectors in our classical solutions, A® —* A® + SA^. 
Expanding the vector Lagrangian in (JTJ) to second order in dA^ and neglecting fluctuations 
in the gravitational field I get, for one of the triad vectors, 

L = -\5F,JF^ - ^A,5A" ~ 2^(^M„) 2 . (12) 

Terms linear in the perturbations vanish if A atl satisfies the classical equations of motion. 
Note that in addition to the mass term proportional to dV/dA 2 , there is an additional con- 
tribution proportional to d 2 V/d 2 A 2 that breaks Lorentz invariance (because of the coupling 
of SA^ to the classical, non-vanishing vector A^.) Therefore, the quantization of SA^ is 
expected to be significantly different from the one of the tachyonic vector in the Lagrangian 



III. COSMIC EVOLUTION 



Our next task consists in studying the evolution of the cosmic triad in a universe that 
contains additional forms of "matter", like an inflaton 4 (pi n f ~ — pi n f), radiation (p r = p r /3) 
or dust (pa = 0). Ideally, we would like the cosmic triad to remain subdominant during most 
of cosmic history, and just around redshift z s; 1 come to dominate the energy density of 
the universe and trigger a stage of accelerated expansion. 

The vector equation of motion (jUJ) is formally the same as the one for a self-interacting 
conformally coupled scalar field. Indeed, the term in parenthesis is proportional to the Ricci 
scalar R, which vanishes during radiation domination. Of course, the similarity between the 
equations of motion of a vector and a conformally coupled scalar arises from the conformal 
invariance of the Maxwell Lagrangian. In some instances it is going to be more convenient 
to deal with a set of two first order differential equations, rather that with a single second 
order one. Introducing the number of e-folds dN = dloga as a time variable, and defining 

B = A + HA, (13) 

the vector equation of motion © can be recast as 

dB ^ n 1 dV n MA , 

dN +W+ HdA=°> (14a) 

^-Z +A = , (14b) 
dN H K J 



where the Hubble constant is given by equation (|7a|). 

In the following I study the vector equations of motion in two limits: the limit where 
matter dominates (early stages of cosmic evolution) and the limit of triad domination (late 
stages). The evolution of the triad depends on the self interaction V. In this paper I focus 
on a class of run-away power-law potentials 

/ A 2 \ ~ n 

V (A 2 ) = M 4 ■ (—) , (15) 

where M is a positive parameter with dimensions of mass. This class turns to be sufficiently 
general for our purposes. In five spacetime dimensions, a self-interaction of this form (with 
n = 1/2) leads to inflation along "our" three spatial dimensions, while keeping the size of 



The inflaton is the component or components responsible for an eventual stage of inflation. 
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3511 . These runaway potentials are also 
lo j. Note that for n > 0, the case we 



the remaining fifth dimension essentially constant 
reminiscent of a class of "tracker" quintessence 1C, 
are interested in, all these models are tachyonic. Non-tachyonic interactions do not appear 
to be particularly interesting. 

A. Matter domination 

Suppose that the energy density of the universe is dominated by the inflaton, radiation 
or dust, i.e. the contribution of the triad to the total energy density of the universe is 
negligible. The scale factor then grows like 

tV , 2 



a (t) = a o — , where (3 = — — (16) 

VW 3(1 + w m ) 

and w m is the matter equation of state. Thus, /3 > 1 during nearly de Sitter inflation, 
(3 = 1/2 during radiation domination and (3 = 2/3 during dust domination. We shall keep 
(3 as a free parameter and consider solutions of the equation of motion where A also grows 
as a power in cosmic time. The reader can verify that indeed, 

A = ( Mn + I) 2 \ ^ (J_\ *k {m 

\3(3(n + l)-n + (3(2(3-l)(n + l) 2 J \M~ l ) 1 ; 

is a solution of the equation of motion (J^J) for an interaction given by equation (|15p. Actually, 
this solution is also an attractor. Perturbing A — > A + 5 A and linearizing equation © for 
the given unperturbed background ([To*]) I find 

/ a\ d 2 V 
5A + 3H5A + (H 2 + -)5A+—5A = 0. (18) 

The term in parenthesis is positive during dust domination, and it vanishes during radiation 
domination. At the same time, the derivative d 2 V/dA 2 is always positive for the potentials 
()15|) . Therefore, small deviations from the solution (JT7j) oscillate and decay away. Note that 
along the attractor, the equation of state of the triad is constant, 

ii 

1 +w A = —— (l + w m ). (19) 
1 + n 

ience, these solutions are analogous to the quintessence and trackers discussed in [l^ and 
lo| . The triad equation of state is always smaller than the one of the dominant component, 
and the closer the latter to —1, the lesser is their difference. The triad equation of state as a 



n 



w, 



0.5 



1 



1 . 5 



2 



'A 



-0 . 6 



-0 . 4 



0.8 



0.2 



-1 




FIG. 1: A plot of WAi equation (|19|). for (from top to bottom) (5 = 1/2 (radiation domination), 
(3 = 2/2) (dust domination) and /? = 15 (nearly de Sitter inflation). 

function of n is plotted for various values of (3 in figure [TJ Because at the present epoch the 
universe is not dominated by dust any more, the value of the triad equation of state today 
lies between the ones along the matter attractor and the de Sitter attractor I will discuss 
next. Hence, experimental constraints 

HQ 

on the value of wde restrict the possible values 
of n. In particular, in order to obtain wa ~ 0.75 today, models with n ;$ 1/2 have to be 
considered. This restriction applies provided initial conditions are chosen in the basin of 
attraction of the tracking solution. 

For later purposes, let me also discuss an approximate solution of the system (JHJ). Sup- 
pose that B <C HA and dV/dA <C HB. Loosely speaking, these inequalities are attained in 
the limits of large A or large H. Then, equations (JTlJl have the approximate solution 



Along this solution the kinetic energy of the field B 2 oc e~ 4N decreases, whereas the potential 
energy V cx e 2nN increases. Hence, soon the potential energy dominates the kinetic one, 
B 2 <C V, so that along this approximate solution the equation of state is 



which is less than —1. However, this approximate solution only holds temporarily, as the 
assumptions we have made finally become violated. Later on I will discuss an example where 
initial conditions place the triad on the approximate solution (J2T1J) for a significant period 



A = A {0) exp[-(iV - N )}, B = B {0) exp[-2(iV - N )) 



(20) 



wa ~ — 1 n, 



(21) 
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of time. Note that in the opposite limit, the limit where the kinetic energy dominates the 
potential one, B 2 ^> V, the triad equation of state is radiation-like, wa ~ 1/3. In both 
limits, the triad is quite different from a (canonical) scalar field (with positive potential), 
for which — 1 < w < 1. 



B. Triad domination 

Along the attractor (J17)) . the energy density of the triad decays slower than the one of 
matter. Therefore, sooner or later the triad will come to dominate the universe. Consider 
thus the equation of motion (jUJ) when the triad dominates the content of the universe. Using 
equations (0) and Q the vector equation of motion © takes the form 

■x orr; dw n dW 4ttG A nT , JV \ . dV 

A + 3HA + lA= ' where 7£T S — { 12V - 3 dA A ) A + lI' < 22 > 

which is just like the equation of motion of a minimally coupled scalar with appropriate 

potential (we won't need the explicit form of W(A)). 

Solutions of equation ()22j) with constant A = A* can be easily identified by requiring 

dW/dA to have a zero at A*. This leads to the condition 

d\ogV 



d\og(4nGA 2 - 1) " (2 '"' 
Those constant A solutions are expected to be stable (attractors) if d 2 W/dA 2 > 0, which 
implies 



dV 1 d 2 V 



+ : ^ ,n t J l-^GA 2 ) 



dA 2 47rGd 2 A 2 



> 0. (24) 



A, 

Along the attractor the energy density is, from equations (ffaj) and (jHaj) . 

W{Al) 

PA = 1-^GAV (25) 
i.e. a constant. Hence, these constant A solutions are de Sitter attractors with wa = — 1; 

they are natural candidates to accommodate the present accelerated expansion of the uni- 
verse. Figure El shows a phase diagram of the system (}2*2*|) when the triad is the dominant 
component of the universe. Note that along the de Sitter attractor, the kinetic energy of 
the triad does not vanish. 

Inserting the potentials (|15|) into equation (J23)) and verifying condition (J24j) I find that 
there is a single (stable) de Sitter attractor at 

Al = J: (26) 
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FIG. 2: A phase diagram of the system (|14|) for n = 1/2. The units are arbitrary. The diagram 
clearly shows that all phase trajectories converge to the single de Sitter attr actor at constant A 
and H. 

Substituting this value of A* into (|2*5|) and requiring it to be of the order of the present 
energy density, one can then estimate what is the required value of M to fit the present 
stage of accelerated expansion, 

M ~ 10"^G- 1/2 . (27) 

Therefore, for n ~ the required value leads to the "infamous" scale M ~ 10 -3 eV, whereas 
larger values of n result into more reasonable energies. Note that we are fitting, rather 
that explaining, the time cosmic acceleration begins. As mentioned above, the exponent n 
determines the value of the equation of state today. The close n to 0, the closer wa is to —1 
today. 




C. Two Examples 

In this section I present two particular examples of possible realizations of vector dark 
energy. Current limits on the value of the equation of state of dark energy wde and its 
derivative with respect to redshift dwoE/dz at z = are 

^U = -1.31±°;i and ^|^ = 1.48±8:tS, (28) 

where the prior on the density parameter of dust VLm = 0.27 ± 0.04 has been assumed 
Because the uncertainties are correlated, the reader is advised to look at the constraints on 
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FIG. 3: Triad energy density (in units of today's critical energy density) versus log 10 (l + z) for the 
interaction (|29|) . Shown is the energy density of vector dark energy for four sets of initial conditions 
(continuous lines). For reference, the energy densities of radiation (dash-dotted) and dust (dotted) 
are also displayed. Note that despite the big difference in the initial value of the energy density (50 
orders of magnitude), the present value of the dark energy density parameter still is Qde ~ 0.72. 



the wde — dwDE/dz plane in figure 10 of |l|. Note that the limits (j2Hj) are significantly 
weaker than the ones derived assuming that Wde is constant. 
The first example has a potential 



and I am assuming that the Hubble constant today is H Q = 70 km/s/Mpc. Note that in this 
model V is not an analytic function of A 2 . In particular, A 2 is constrained to be spacelike. 
Figure |3] shows the evolution of the energy density for different sets of initial conditions, 
and figure 0] shows the corresponding behavior of the triad equation of state. As clearly 
seen in the figures, the late time evolution of the cosmic triad is quite insensitive to initial 
conditions. For this particular model, Wa ~ —0.87, dwA/dz « 0.10 and Qa ~ 0.72 today. 
These values are consistent with the current limits on the properties of dark energy ([28)1 . see 
also figure 10 in [1]. Larger values of n lead to violations of the limit on wde today, whereas 
smaller values yield was closer to —1. 

In the second example the interaction term is given by 




where 




(29) 




1/2 




3/2 



where 



(30) 
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FIG. 4: A plot of the triad equation of state for the same initial conditions and interaction as in 
figure |31 The attractors are easily identified by their constant value of wa, which can be read off 
figure n For ° ne set of initial conditions, the system reaches first the radiation attractor, then 
proceeds to the dust attractor and finally to the de Sitter attractor. For another, the system does 
not reach the radiation attractor, but does reach the dust one before continuing to the de Sitter 
solution. And lastly, for yet a different one the universe barely reaches the transition to the de 
Sitter attractor just in time. For all initial conditions, wa ~ —0.87 today. 



This corresponds to n — 1/2 in equation (|15p. In this case the triad evolution can be 
consistent with current observations if the initial value of A is fine-tuned. The initial value 
of B = (A + HA) can be quite arbitrary though. Cosmic evolution for such a model and 
tuned initial conditions is shown in figures 03 and El The different lines correspond to different 
initial values of B. Note that all of them yield the same final result, and in particular all 
of them join the approximate solution (|2*U|) . As seen in figure El the equation of state along 
that solution remains constant at Wa = —4/3 (equation ()21Jl) all the way until today. By 
construction, in this particular example wa ~ —1.33, dwA/dz ~ and Qa = 0.72 today. 
For integer values of n and fine-tuned initial conditions it is also possible to obtain histories 
where the equation of state wa remains constant at a value significantly smaller than —1 for 
a long period of time, and just recently approaches —1. For these models dwA/dz tends to 
violate the limit in (}2*B|) . Certain analyses of supernova data ^| suggest that the equation 



of state has evolved from Wde ~ at z 1 to Wde < — 1 today 



14j. Within the class of 



models (|15|) I have not been able to obtain such a behavior while keeping at the same time 



14 



50 - 



\ 40 




FIG. 5: Triad energy density (in units of today's critical energy density) versus log 10 (l + z) for 
the second model (|5U|) . For reference, the energy densities of radiation (dash-dotted) and dust 
(dotted) are also displayed. The shown trajectories share the same initial value of A, though 
the initial values of B are different. The present value of the dark energy density parameter is 
tt DE ~ 0.72. 
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FIG. 6: A plot of the triad equation of state of the dark triad for different initial conditions (same 
as in Figure I2J). For all of them the triad reaches the equation of state wa = —4/3, equation (|21|) . 
and remains at that value till the present epoch. 

an early epoch of radiation domination. 



15 



IV. PERTURBATIONS 



The properties of any dark energy candidate not only comprise its equation of state, but 
also the way their perturbations (if any) behave 8}. These perturbations are coupled through 
Einstein's equations to metric perturbations, which in turn affect observables like CMB 
temperature anisotropies. The evolution of triad perturbations can also help to determine 
whether models with wa < —1 suffer from serious instabilities and are hence unviable. 
Unfortunately, cosmological perturbation theory with a triad turns to be rather involved, 
as scalar, vector and tensor modes couple to each other. Therefore, in this section I only 
scratch the surface and mainly present qualitative features of the equations. 

In scalar longitudinal gauge and vector gauge, the most general linearly perturbed spa- 



tially flat FRW metric is 
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42, 



43] 



ds 2 = a 2 (7]) {-(1 + 2$)dr/ 2 - 2B X d^dx* + [(1 - 2*) 6 {j + 2h ij ]dx i dx j } , (31) 

where B % is a transverse vector, diB % = 0, and hy is a transverse and traceless tensor, 
dih % j = h l i = 0. In a similar way, we can decompose the perturbations of the vector fields 
in the triad 5A^ into scalar and vector components, 

5Al = (8Ald iX a + 5A«). (32) 

Here 5Aq and \ a are scalars and SAf is a transverse vector, di5A m = 0. Note that for 
convenience we are using conformal time. In the following, spatial indices are raised and 
lowered with the metric 5ij. 

A. Linearized triad equations 

Because the vector fields in the triad do not couple to each other, we can study each of 
them at a time. Let me denote by Af the spatial components of one of the background triad 
vectors and let me drop in the following the index a. The 0-component of the linearized 
equation of motion (jSJ) in the spacetime (j3Tj) is 

- 2a 2 %) 6A = —BiA 1 " - + V)A 1 ' + d^x', (33) 



dA\ 

where a prime denotes a derivative with respect to conformal time. Remarkably, even 
though this is a scalar equation, it contains the vector perturbation B l . In that respect, 

16 



non- vanishing vector fields lead to violations of the decomposition theorem 17[. In the 
absence of a background vector quantity, the only way to obtain a scalar linear in a vector 
B i is to compute its divergence d^B 1 , which vanishes for a transverse vector. However, if 
the spacetime contains a non-vanishing background vector it is possible to construct a 
scalar linear in a vector perturbation through the combination AiB % . The reason why the 
decomposition theorem is expected to hold in FRW universes is that in general they do not 
have any spatial direction singled out. 

The i-components of the linearized equations have the form of the gradient of an equation 
that involves scalars plus an equation that involves transverse vector quantities. The scalar 
equation is 

X " -6A' -6 = -2a 2 ^ X , (34) 
where 5 is implicitly defined by the decomposition 

2®A'! + ($ + - 2a 2 ^S(A 2 )A t = + 6 t , (35) 
into a scalar 5 and a transverse vector 5j. The perturbation of A 2 is given by 

6(A 2 ) = 2a- 2 [£{$6^ - h l3 )AJ + A\d iX + SAA] , (36) 

and in particular, does not contain 5A . The remaining vector equation takes the form 

6 A! - djd j 5Ai - (A'd^Bi - 2h^' - ^ = -2a 2 ^-5Ai. (37) 

Therefore, again, scalar, vector and tensor perturbations are not decoupled, since there 
exists a non-vanishing vector in the background. 

Note that ()33|) is not a dynamical equation for d~A , but a constraint. It seems that for 
tachyonic models, dV/dA 2 < 0, one cannot solve for 5Aq, since in that case the operator 
did 1 — 2a 2 dV/dA 2 is not invertible. This would point out to a potential inconsistency of 
tachyonic models. However, as we shall see next, this difficulty is only apparent. Taking the 
Laplacian of equation (jMj) and using the time derivative of the constraint (}3*3"|) to substitute 
the value of did l 5A' , one obtains a first order differential equation in time for SA that does 
not contain its spatial derivatives, 



dA 5A ° 



dV\ a' dV 
+ 2- 



dA 2 J a dA 2 



5A = f($,*, X ,a,A,B). (38) 
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Here, / is a function of the specified variables and its derivatives we shall not be concerned 
with. What is important is that if a set of initial conditions satisfying the constraint (}3"3"j) 
is specified, then, the solution of equations (J3*8|) and ()34|) is guaranteed to satisfy equation 
(1331) at all times. Hence, it is in fact possible to solve for SA . 

To conclude this part let us count how many "degrees of freedom" (per vector) the 
triad perturbations contain. We have seen that 8Aq is constrained, so it is not dynamical. 
Equation (J3"4"|) contains the second time derivative of % (1 dof) and equation (|37j) the second 
time derivative of the transverse vector 5 A (2 dof). Thus, there are 3 degrees in freedom in 
total, as pertains to a massive vector. 



B. Perturbed Energy Momentum Tensor 

The perturbations in the triad induce perturbations in its energy momentum tensor, 
which in turn are responsible for sourcing metric perturbations. In this subsection I shall deal 
with vector and tensor perturbations, which cannot be sourced in conventional cosmological 
models. 

The perturbations in the energy-momentum tensor of the triad vectors A ^ can be de- 
composed into an isotropic pressure 5p A and a traceless anisotropic stress IPj perturbation, 
ST l j = Sp^j+pyJ^j- The anisotropic stress itself can be decomposed into scalar, vector and 
tensor components, IP, = ^Wj + ^Wj + (l,) n/ + WlPj-, where WlT, = (9*9,- - 5 i j d k d k /3)U, 
Wff; = c^n* (TP transverse), and is transverse and traceless. 

fin 

The equations of motion for vector (metric) perturbations are p, LL2( 

djB*' + 2-d j B i = STrGaW^n'j, (39) 

where ( v 'H i j is the vector part of the triad anisotropic stress tensor. In the absence of 
anisotropic stress sources, equation (J3"9"|) implies that vector perturbations decay away, 
B l oc 1/a 2 . So even if they are generated in the early stages of the universe, they are 
not expected to be significant today. The transverse and traceless anisotropic stress sources 
tensor perturbations, i.e. gravitational waves, 

K" + 2^ V - d k d k h ij = 8nG PA {t) U ir (40) 

As opposed to vector perturbations, in the absence of sources the amplitude of long- 
wavelength gravitational waves remains constant. Hence, if they are primordially produced, 
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say during an inflationary stage, they could still have a sizable amplitude today. 

I shall not write down all the terms that the perturbed spatial components of the energy 
momentum tensor of a single triad vector contains. They straightforwardly (but tediously) 
follow from the insertion of the perturbations (|31|) and (|32|) into equation (J2J). Instead, for 
the sake of illustration I shall consider only 

5["r*j] = 2a- 2 ^- 2 (A'SAj + 6A%) + • ■ • , (41) 

where the dots denote the multiple terms I am not explicitly writing down. In order to 
study the evolution of vector and tensor perturbations we have to compute the vector and 
tensor components of the previous expression. In Fourier space, these are given by 

PA Wn* = -^^!^(fc.^)M- (42) 



a 



dA 2 

a 



ui k \ k ■ A a 

A m 5A a , + 8A ai A a , - ( —6 A* + 5A ai ^ — — 

x k J k I k 



- ' 5 ) ~ I A a ■ 5A a 



(43) 



which as required are transverse and traceless. Note that we sum over the three triad vectors 
to obtain the total energy momentum tensor perturbation. Because the triad perturbations 
SA ai are a priori totally independent from each other, the anisotropic stress WD?j does not 
vanish in general. Thus, vector fields are not only expected to source vector perturbations, 
but also tensor perturbations. Again, the reason is that the decomposition theorem is 
violated. With the aid of the background vectors it is possible do construct traceless and 
transverse quantities linear in the perturbations. 

If the perturbations have certain symmetries , does indeed vanish (for the particular 
term in the energy momentum tensor we are considering). Because A a% and k % are the only 
vectors in the problem, the triad perturbation 5A ai is expected to be a function of A ai and 
k\ Since 5A a is transverse, it has to be of the form 

6A ai = a ■ (Vj -jj^j A aj + (3 ■ e^kA^k*, (44) 

where a and (3 are two scalar functions and eyjt is totally antisymmetric. Assume that a and 
(3 do not depend on the index a (that is, they do not depend on kiA ai .) Then, substituting 
equation in (|4^j) and using (0) one finds not only that = 0, but also ^IT = 0. 
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C. Stability 



The main worry one faces when dealing with tachyonic fields is their stability. In order 
to figure out to what extent the solutions we have found in Section ITTT1 are stable, we should 
solve the system of cosmological perturbation equations just presented. Obviously, the 
coupling between scalar, vector and tensor perturbations makes this task quite formidable. 
In this section we dramatically simplify the equations by neglecting metric perturbations and 
concentrating on a particular set of vector modes. The hope is that this drastic simplification 
captures the qualitative features of the equations. 

So let's set $ = ^ = B l — hij = and consider the triad perturbation equations in 
Fourier space. From equation (J35j) . for modes for which k ■ A — 0, it follows that 5 = 0. 
Then, x — ^A = is a solution of equations (}3"3"j) and (J3"4~j) . Because 5i oc (A^5Aj) Ai, if in 
addition we consider perturbations such that 5 A ■ A = 0, the remaining vector equation (jSTj) 
reads 

k 2 dV 

5Ai + HSAi + —5A { + 2—dAi = 0, (45) 
a dJ\ 

where for convenience I have gone back to cosmic time. Note how 2dV/ dA 2 plays the role of a 
mass term in the last equation. This is why I call interactions with dV/ dA 2 < "tachyonic" . 
Generically, if dV/dA 2 is negative, one expects growing modes, that is, instabilities. 

In order to check the stability of the triad, it suffices to consider long-wavelength modes, 
k = 0, in equation ()45|) . For sufficiently high k, the gradient dominates the interaction and 
solutions are stable. Hence, any form of instability is an "infrared" effect, rather than an 
ultraviolet one. For a given expansion, equation (fTHjl . along the attractor (fTTjl . the long- 
wavelength solution of equation (|43|) is 

J\A. i n 

-^ = C + 1 7 + + C_ t 7 - , where 7+ = and 7 _ = -3/3 + . (46) 

a n + 1 n + 1 

Here, C + and C_ are two integration constants and for convenience I have divided by the 
scale factor to obtain the length of the perturbation. We should compare these solutions 
with the behavior of the background, equation (fTTj) . Recall than A is the common length 
of the background triad vectors. The C + mode in equation (|46J) grows as fast as A and the 
C_ mode grows less rapidly than A if n(l — 3/5) < 1 + 3/3. Therefore, the C_ mode decays 
relative to A during inflation, radiation and dust domination. Hence, within the scope of our 
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analysis, the system is not unstable throughout the period where the triad is sub dominant 5 . 
Note however that the triad would be unstable for certain values of n if there had been a 
period of cosmic history during which (5 < 1/3. 

When the triad becomes dominant, the vector evolves towards the de Sitter attractor 
(J2SJ). The solution of equation along the de Sitter attractor is 

XA. 

— = C + exp(Hj) + C_ exp(-2HJ), (47) 
a 

where if* is the value of the Hubble constant along the de Sitter solution. Note that 
along the latter A itself is constant. Thus, the de Sitter attractor is unstable, in the sense 
that 5Ai/a (the length of the vector perturbation) grows relative to A (the length of the 
background vector). Because during the previous stages of cosmic history (when the triad 
was subdominant) the relative amplitude of the perturbations 5 A/ A has remained unaltered, 
the time perturbations become relevant will depend on early universe initial conditions. If 
the primordial vector amplitude agrees for instance with the (scalar) amplitude of density 
fluctuations, 5A/A « 10~ 5 , 5A/A becomes of order one about 12 e-folds after the onset of 
triad domination. By then the universe has presumably grown anisotropic, because there 
is no reason to expect that perturbations in the three different triad vectors are correlated. 
We'll have to wait several billion years till that happens. At present the effect is still small, 
since the relative amplitude has increased at most by a factor of order one. In fact, it is 
tempting to speculate whether some of the anomalies observed in the CMB radiation 
specially in the quadrupole and octopole (see however [44]) and the related hints of statistical 
anisotropy [22], might be due to such an instability, which sets in once the universe starts 
to accelerate and mostly affects large scales. 



V. SUMMARY AND CONCLUSIONS 



In this paper I have considered whether a vector field could be responsible for the current 
stage of cosmic accelerated expansion. The existence of a vector with non-vanishing spatial 
components turns to be compatible with the isotropy of a Friedmann-Robertson- Walker 
universe provided the vector is part of a "cosmic triad" , i.e. a set of three identical vector 

5 Strictly speaking it is not stable either, as the relative perturbations do not decay. 
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fields pointing in mutually orthogonal directions. A set of three identical self-interacting 
vectors naturally arises for instance in a gauge theory with 577(2) or 5*0(3) gauge group. 

A distinctive property of a cosmic triad is that its equation of state of can become less 
than —1, even though its kinetic terms have the conventional form. The necessary condition 
is that the self-interaction is "tachyonic", i.e. it naively gives rise to a negative squared 
vector mass. Although the simple analogy with scalar tachyons suggests that the study of 
tachyonic vectors is justified, there are also arguments that connect tachyonic vectors to 
phantom particles 3s|. In analogy to tracking quintessence models ^j], in this paper I have 
explored a particular class of tachyonic models, where the interaction is an inverse power- 
law. For appropriate initial conditions, there exist seemingly viable cosmologies where dark 
energy has an equation of state w < — 1 throughout cosmic history, including today. In top 
of that, these models have attractors that render cosmic evolution quite insensitive to initial 
conditions and yield equations of state sufficiently close to —1 today. The experimental 
constraints on the value of the dark energy equation of state restrict the power of the vector 
self-interaction, while the time at which cosmic acceleration sets in determines its energy 
scale. Within the class of models I have considered, the current limits on the equation of 
state of dark energy require the vector self-interaction to be a non-analytic function of its 
squared length. 

Finally, I have scratched the surface of cosmological perturbation theory in the presence of 
a cosmic triad. The most remarkable feature is the violation of the decomposition theorem. 
In the presence of a vector which has non-vanishing spatial components, scalar, vector 
and tensor perturbations do not decouple from each other. In particular, metric vector 
perturbations show up in scalar equations, and in the absence of particular symmetries in 
the perturbations vectors are able to source tensors. The (scalar) time components of the 
triad vectors are not dynamical, i.e. they are constrained. Despite an apparent difficulty, it 
is possible to solve the constraint also for tachyonic models. I have also considered solutions 
of the perturbation equations for a particular set of modes, under the assumption that metric 
perturbations are negligible. During inflation, radiation and dust domination the relative 
perturbations in the triad remain constant. However, there is a long wavelength instability 
during the late-time stage of de Sitter acceleration, where triad perturbations grow relative 
to the background. Hence, the time the universe becomes anisotropic depends on early 
universe initial conditions. For reasonable primordial perturbation amplitudes, the universe 
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is expected to become anisotropic long time after the onset of cosmic acceleration. The 
instability of the triad during this epoch also suggests a possible relation between the large 
angle anomalies in the CMB sky and the onset of cosmic acceleration, but further work is 
needed to test this idea. A more careful investigation is also needed to establish whether 
tachyonic vectors are fully stable, and how the inclusion of metric perturbations affects the 
behavior of the triad perturbations (and vice versa). 

To conclude, at the level of the present analysis it seems that vector fields could indeed 
be responsible for the present stage of late time cosmic acceleration, though it is yet unclear 
how quantum mechanics constrains the tachyonic models I have studied here. 
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APPENDIX: A NON-ABELIAN TRIAD REALIZATION 

In this appendix, I show that a cosmic triad can arise from a non-Abelian SU (2) gauge 
theory. For that purpose, consider the gauge-invariant Yang-Mills Lagrangian 

1 



-- J d*xy/=gF^Fjr, (A.l) 

where F^ v is the non-Abelian field-strength 

F; v = d,Al - d v Al + s\ c A\Al. (A.2) 

The totally antisymmetric tensor e a b c encodes the structure constants of the SU(2) Lie- 
algebra. The equations of motion of the fields in an arbitrary spacetime are 
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=d„ {V=gF a ^) + e\ c A\F^ v = 0. (A.3) 

The ansatz (jSJ) satisfies these equations of motion in the FRW universe (HJ) if 

A + 3HA + ( H 2 + - J A + = 0, (A.4) 
\ a J dA 
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where 

V = \ (A,A») 2 • (A.5) 

Equation (|A.4|) agrees with the triad equation of motion The energy momentum tensor 
of the gauge fields is 

T, v = F; p F a S - \F a pa F^g, u . (A.6) 
It can be verified that T°j = 0, whereas 

-T° = |(i + iM) 2 + 3V, = Q(i + #A) 2 - 3V + 2^A 2 ^) 5V (A.7) 

The self- interaction V is again given by equation (|A.5|i . Therefore, the energy density and 
pressure of the non-Abelian gauge fields agree with the ones of the triad, equations 0- This 
equivalence between the triad and the non-Abelian gauge fields in the symmetric case we 
are considering is confirmed by substituting the ansatz (JHJ into the actions Q and (jA.ljl . 
Note that for the triad vectors, as opposed to scalar fields or perfect fluids, the Lagrangian 
density is not the pressure nor the energy density. The existence and some properties of 
non-trivial solutions of general relativity coupled to an SU{2) gauge field in a FRW universe 
have been considered in 1451 . 
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